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We consider a fibrillar medium with a continuous distribution of two-level atoms cou- 
pled to quantized electromagnetic fields. Perturbation theory is developed based on the 
current algebra satisfied by the atomic operators. The one-loop corrections to the disper- 
sion relation for the polar itons and the dielectric constant are computed. Renormalization 
group equations are derived which demonstrate a screening of the two-level splitting at 
higher energies. Our results are compared with known results in the slowly varying enve- 
lope and rotating wave approximations. We also discuss the quantum sine-Gordon theory 
as an approximate theory. 

PACS: 42.50.* 



5/96 



1. Introduction 

The interaction of radiation with a medium of atoms is an important problem with 
many applications, in particular to the propagation of optical pulses and to lasers. When 
studying resonance phenomenon it is natural to approximate the atoms as two-level sys- 
tems, resulting in optical Bloch equations. For a single atom in an external classical electric 
field, this problem was solved by Rabi[0. A collection of atoms coupled to radiation is de- 
scribed by the so-called Maxwell-Bloch equations. Propagation effects in the latter model, 
were studied by McCall and Hahn]^]. There, a semi-classical approximation was made 
wherein the electromagnetic field was purely classical. With some additional approxima- 
tions (see below) it was shown that the equations of motion reduce to the well-known 
sine-Gordon equation, and the optical soliton solutions were observed experimentally also 
in . The problem of many-atom spontaneous emission was studied in a simplified model 
by DickeP], which does involve quantized electric fields. We refer the reader to 0] for an 
excellent account of optical resonance phenomena. 

In this paper we study the fully quantum system of radiation in interaction with a 
continuous distribution of atoms in a fiber geometry. Two different models are studied, 
one following from taking the interaction hamiltonian to be —d ■ E, which we refer to as the 
'current-model', the other following from the minimal coupling prescription p ^ p — eA and 
referred to as the 'charge-model'. We develop perturbation theory using the current algebra 
satisfied by the atomic operators. This allows us to easily determine the dependence of 
various physical quantities on the number of atoms A^. In particular by computing the 
photon self-energy we determine the first quantum corrections to the polariton dispersion 
relation and dielectric constant. 

In the current model, we show how the quantum corrections imply a renormalization 
group equation for the two-level energy splitting, this splitting becoming screened at higher 
energies. 

In the last sections of the paper we compare our results with known results obtained 
under various approximations, namely the slowly-varying envelope and rotating wave ap- 
proximations. We also argue that the quantum sine-Gordon theory has some validity as 
an effective quantum field theory. 

The two-dimensional quantum field theories we study are interesting in their own 
right. The more interesting of the two models we study (the current- model) is defined by 
the hamiltonian 
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where Hq is a free hamiltonian for a scalar field 0, ^'3, satisfy a current algebra (see 
equation 2.29), luq is the two- level splitting of the atoms, and /? is a dimensionless coupling 
that depends on the strength of the dipole transition. The renormalization group equation 
for the radiative level shifts is a consequence of the beta-function for u>o, which we compute 
to lowest order. To the best of our knowledge, this model is not integrable, though various 
approximations to it are integrable (see below). 

2. Two Models for the Quantum MaLxwell-Bloch Theory 

There are two related models describing the coupling of two-level atoms to quantized 
electromagnetic fields. One follows from taking the interaction hamiltonian to be —d ■ E 
where d is the dipole moment operator and E the electric field. The other follows from the 
usual minimal coupling prescription, p-^p — eA where p is the momentum operator and 

— * 

A the vector potential. We will refer to these as the 'current- model' and 'charge-model' 
respectively. In this section we describe the reduction of both these models to one spacial 
dimension. We will first consider the case of a single atom, and then extend this to a 
continuous distribution of atoms. We set c = 1 in most places, but keep ^ 7^ 1 in some 
formulas to clarify certain points. All formulas with 1 have both c and h restored. 

For both cases, we model a single unperturbed atom as a single electron which has 
two eigenstates |2), |1), with energy difference fiujo- Letting Hq^°^ denote the unperturbed 
atomic hamiltonian, one has 
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(2.1) 



In the basis (|2),|1)), 
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2.1 Current-Model 

One can couple the 2-level atoms to the electromagnetic field starting from the energy 
of an electron in an external electric field. Let the single atom be centered at a^o, and let 
d = ex denote the electric dipole moment operator. The complete hamiltonian is 

H = H^''^'^ + ifo^*"'" - d- E{x), (2.3) 

where Hq^^^ is the free Maxwell hamiltonian, and x is the position of the electron. Let us 
further assume that the electric field does not vary significantly over the region where the 
atomic wavefunction is non-zero, i.e. 

{a\d- E{x)\b) ^ E{xo){a\d\b), (2.4) 

for \a,b) unperturbed atomic eigenstates. 

— * 

The matrix elements of d have the following general form 

(l|dll) = (2|dl2) = 

2.5 

{2\d\l) = de"^n, {l\d\2) = de-"^n, 
where d is a real parameter and n is a unit vector that specifies the orientation of the atom 

— ♦ — * 

in space. {a\d\b) — ior a = b since d is a vector operator with odd parity and the states 
I a) are assumed to have definite parity. 
Letting E = E -n, one then has 

d- E{x) = dE{xo) (e*"cT+ + e-'^d") , (2.6) 

where 

a+ 

We will later need the algebra of the u-operators: 

[(73, a±] = ±2(7±, [a+, a-] = (73. (2.7) 

Note that the phase e*" can be removed by letting a"^ ^ e^'^°'a^ without affecting the 
algebra; we henceforth set this phase to one. 

Consider a fibrillar geometry, where the atom can be viewed as an impurity in an 
optical fiber of length L and cross-sectional area where L » \fA. One can perform a 
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reduction to this essentially one-dimensional theory as follows. The free Maxwell action 



which determines Hq^^'^ is 



^ j d^xdt , (2.8) 



'S'Maxwell . 

where Ff^i, = df^A^, — d^A^j,. Let x denote the direction along the fiber, and y,z the 
directions transverse to it. One can first consistently set = 0. We also require the 
energy flux to be along the fiber, so that the Pointing vector is in the x direction. This 
requires A^ = ^ and dyAz — d^Ay — 0. Thus, we only have to deal with the components 
of A transverse to the direction of the fiber. Of these, only A = A-n couples to the atom. 
Assume A is independent oi y, z, and let J dydz — A. One has 

^Maxwell = ^ / ^xdt^ (dtAOtA - d^Ad^^Aj . (2.9) 

One can alternatively understand the appearance of the cross-sectional area A by 
considering the mode expansion of a free scalar field in finite volume V, 

^(^) = E -7t7^ (a{k)e-'^-' + a\k)e'^'-^) . (2.10) 

In the fibrillar geometry, V = AL and VA is very small compared to L. As L — > oo, 
the modes in the ^-direction are nearly continuous. The transverse modes in the y, z 
directions have a maximum wavelength on the order of vC4 and are thus very energetic 
in comparison to the low energy modes in the L-direction. Thus, we are assuming these 
high energy transverse modes are negligibly excited, which is reasonable if l/\fA ^ a^o- 
In general A should be replaced by ^eff, which is the effiective cross-sectional area of the 
fiber as a waveguide. 

Rescale A: 

Then 



'S'Maxwell 



= j dxdt^ idt(j)dt(j) - d^(j)d^(j)) . (2.12) 



n 

The field is a dimensionless scalar field. In the quantum theory it satisfies the commu- 
tation relations 

[(j)(x, t), dt^{x', t)] = i5{x - x'). (2.13) 



Using E = —dfA, and dividing by h to give the hamiltonian units of 1/time, we obtain 
the complete hamiltonian of the current-model: 

H = H^ + H^'""^ + , (2.14) 

where iifg*"™ is defined in 



hS = I dx ^[dtcl^f ^l{d^cl^)\ (2.15) 
iyi.™^^ = f 9,<^(^o)K + a-), (2.16) 



and 



The terminology 'current-model' refers to the fact that the spin operators couple to dtcj) 
which is the spacial component of the conserved topological current t^^dj^cj), where e^'^ is 
the anti-symmetric tensor. 

The parameter 13 is the important dimensionless coupling constant of the model. Since 
d ~ ei?atom, where -Ratom is an atomic dimension, 

~ 2 I — — I ~ -'''atom ^2 j^g") 

Stt ^ \hcj AeS ^ 137 AeS ' 

Thus, generally, /3^/87r is very small. An idealized upper limit would correspond to a chain 
of atoms in a waveguide that is one atom in thickness, so that -Ratom ~ -^efF- This is 
perhaps nearly realizable with a polymer waveguide. In this situation /3^/47r ^ 1/137. We 
will refer to this hypothetical limiting case where the quantum effects are strongest as the 
quantum optical chain. 

The parameter (3'^ /Sir determines the spontaneous decay rate 1/r of a single excited 
atom. First order perturbation theory gives 

\{l,k\2)\' = ^{2nS{uo-\k\)f (2.19) 

where |1, k) = |l)atom ® \k) and \k) is a one-photon state with wave-vector k. This leads 
to ^ 

^ = J dk\{l,k\2)\^ i2iT5{uo-\k\))-^ = (^uo. (2.20) 
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2.2 Charge-Model 

For the charge- mo del we begin with the standard way to couple particles to electro- 
magnetic fields and consider a single electron hamiltonian 

H-^^{p- (^MS)f + V{x), (2.21) 

where p is the momentum operator and V is the atomic potential. Again, we make a 
2-level approximation for ifg*"™ = p^/2me + V , and consider states |1, 2) as in (|2.1| ). As 
for the current-model, we assume the analog of (|2.4| ). Using 

Tfl 

p = ~z^[x,H^'°^], (2.22) 



one has 



A(xo) ■ {a\p\b) = UEb - Ea)A{xo) ■ (a| d\b), (2.23) 
rUp n 



where Ea is the energy of the state |a), and d is again the electric dipole operator d = ex. 
Using the parameterization (2.5) for the cT matrix elements, and rescaling A as in 



( p.ll|) , one obtains a hamiltonian of the form (|2.14D , where now 



hS = (/ dx ^{dt<pr + l{d^<f>r^ + <p\xo), (2.24) 

and the interaction term is: 

^^(cha.ge) ^ _,^^^(^^) _ . (2.25) 

We refer to this as the 'charge-model' since (j){oo) — (/>(— oo) is the charge associated to the 
topological current e^'^d^cj). 

The spontaneous decay rate to lowest order for the charge-model is the same as for 
the current-model ( |2.2UD . 



2.3 Continuous Distribution of Atoms 

Consider now a collection of atoms, with the N atoms positioned aX x = Xi, i = 1, N . 
Also, let di = e(x — Xi), |1, 2)i and a{i) denote the dipole moment, 2-level states, and Pauli 
matrix operators for the i-th atom. Since i{l\di\l)i = i{2\di\l)i = 0, and i(l|2)i = 0, the 
matrix elements of di don't depend on Xi. However, in general the atoms have variable 
orientations in space and i(2|(ij|l)i = de'^'^rii, where rii can vary from atom to atom. To 
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simplify the situation, we assume that aU the atoms are somehow ahgned, for example by 
some external electric field or by being embedded in a crystal. Then, 

i{2\d,\l), = de'^^n, i{l\d,\2)i = d e-'^^n. (2.26) 

For the current-model, the interaction becomes 

^^(current) = ^J2^tHx^) (a+ (z) + a" (z)) . (2.27) 

i=l 

Introduce space-time dependent spin operators as follows: 

N 

^«(x) = J2 (^"{i)S{x - X,). (2.28) 

1=1 

These operators satisfy a current algebra 

[S3{x,t),S^{x',t)] = ±2S^{x,t)5{x-x') 
[S~^{x, t), S~{x' , t)~\ = Ss{x, t)5{x — x'). 
The hamiltonian for both the current and charge models takes the form (|2.14| ), where now 



(2.29) 



^atom = ^Jdx Ss{x), (2.30) 
^(current) = ^ f dt<i>{x,t) {S+{x,t) + ^-(x,t)) 

,nr (2.31) 

^^(charge) ^ _.UJoP / ^^^^^^ {S+ix,t) - ^-(x,^)) . 



For the charge-model, the additional term in (|2.24| ) leads to a mass term for the scalar 
field: 

at = \dx Q(a,</>)2 + \{d^<\>f + ^ </>2(x)) , (2.32) 

with 

A^ = 4vr^y^- (2-33) 

In working with the above formulation, one must impose a further condition that there 
is a single electron bound to each atom. For the one-atom operators, this is manifest in the 
two-dimensional representation of the a's, which has the additional relations (a"^(z))^ = 0, 



(0-3(2))^ = 1. Note from the definition (]2]28D that 

(^±(a;,t))V0. (2.34) 
These issues are more easily resolved in a fermionic description, which we turn to next. 



2-4 Fermionic Description 

Consider first the one-atom case. In a second quantized description, one introduces a 
fermion wavefunction It/j): 

= + 62I2), (^1 = + {2\bl (2.35) 
where the 6's are fermion operators satisfying 

{61, 61} = {62, 4} = 1- (2-36) 

The operator b\ (63) creates an electron in the lowest (highest) level. The a operators then 
have the following well-known representation: 

a+ = blbi, a- = b\b2, as = b^ - b\bi, (2.37) 

and the algebra ( p77| ) is a consequence of ( |2.36|) . Since each atom has a single electron, one 
must impose the constraint 

4^2 + b\bi = 1. (2.38) 

This additional algebraic relation leads to the 2-level relations (o"^)^ = 0, (crs)"^ = 1. 
In the multi-atom case, one can define 



S+{x) = bl{x)bi{x), S-{x) = b\{x)b2{x) 
Ssix) = bl{x)b2ix) - b\{x)bi{x). 



(2.39) 



The algebra (|2.29D is then a consequence of the anti-commutation relations: 

{biix), b\{x')} = {b2{x),bl{x')} = 6{x - x'). (2.40) 
Define the number operator: 

N = j dx (bl{x)b2{x) + 6l(x)6i(x)) . (2.41) 

The operator commutes with the hamiltonian and corresponds to the number of atoms 
in the sample. We therefore impose the constraint 

N = N. (2.42) 



In evaluating quantum transition amplitudes, one doesn't have to take special account of 
this constraint as long as one deals with initial and final states that satisfy the constraint, 
since [N,H] = 0. Note that 

J dx Ssix) = -N + 2 J dx bl{x)b2{x). (2.43) 

In the sequel, we will work with the perturbative vacuum: 

|0) = |0)photon ® I Datom, (2.44) 

where | |) denotes the atomic state with all N atoms in their lowest energy state, and |0) 
is the state with no photons. This vacuum state satisfies 

b2{x)\n) = S-{x)\n) = 0, (2.45) 



and 



From ( |2.43| ), one has 



= (2.46) 



j dx Ssix) \n) = -N\fl). (2.47) 



Assuming spacial translation invariance of |0) (we are ignoring boundary effects for L very 
large), ( |2.47|) implies 

s,{x)\n) = -^\n). (2.48) 

We remark that \Q) is not the exact ground state of the theory; see section 5. 



3. Perturbation Theory 

In this section, we study the perturbative expansion for the general correlation func- 
tions of the field (p, which are simply related to electric field correlators via ( p.ll|) . For 
convenience of notation, in the charge-model we make the redefinition i^^S^, which 

does not affect the commutation relations ( |2.29| ). The interaction hamiltonian for both 
models can then be written as 



Hint=g j dx 0{x,t) {S+{x,t) + S~{x,t)) 
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(3.1) 



where for 



current — model 



9=^, = dtc 



charge — model : g 



(3.2) 



O 



We begin with the partition function, defined as the vacuum to vacuum transition 
amplitude. In the interaction picture, 



z = {n\u{oo,-oo)\n), 



(3.3) 



where 

U{oo,-oo) =Texp(^-ig j dxdt 0{x,t){S+ + S')^ . (3.4) 
The correlator (| \S°''^ ■ ■ ■ S"'" \ l)g=o is only non-zero when J2i o-i = 0? thus, 

^ = 5Z(^"^'" / dX2n---dXi / dt2n---dti {0\O{t2n, X2n) ' ' ' O{xi,ti)\0) g=o 

n=0 "J t2n>--->t\ 

X 5^ (i \S'''-{x2nMu)---S''^{xiM) \ i)g=0. 

(3.5) 

The time dependence of the S-correlators is simple, since the hamiltonian at = 
implies 

^±(x,t) =e±*'"o*^±(a;,0). (3.6) 

Thus, 



a |^"^"(a:2n,t2n)---^"H^l,^l)l i)<?=0 



■ 2n 

n 



i=l 



(i |^"^"(a;2n,0)---^«^(xi,0)| i),=o. 

(3.7) 

It will be helpful to pass to momentum space. Our conventions are k = {uj,k), 
d'^k = divdk, k"^ = uj"^ — k"^ , and k ■ x = ut — kx. Define 



2n 5 X2n J 



d^kl d^k 



2n 



(27r)2 {2ny 



' 2n 

n 

\i=l 



) G'^^^)^kuk2, ,M 

(3.8) 



Substituting this into ( |3.5|) , one can perform the time integrals using 



dt'e-'^^' 



-icjt 



uj + ie 



(3.9) 
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where here and below e is infinitesimaUy smaU and positive. One finds 

2n 



X /„,.....„,.. (i^i,...,i^2„-i) Yl W |S''="(i2„,0).--S'"(ii,0)| Ds^o, 



ai,..,a2„=± 



(3.10) 

with 



2n-l 



/ai,...,a2„(Wl,...,W2n-l) = TT — ^ ; ^ (3.11) 

The infinite volume singularities in the sum ( |3.1C1| ) must be regulated in order to obtain 
something meaningful. However, the correlation functions, which are obtained by dividing 
by Z, are well defined order by order in perturbation theory. 

Let $1, $2, ••• denote some local fields. The same analysis as above leads to 

72 U ^2 U 



(o|T($i(x;)---$^(x^))|o) = G(-Hk[,g, ,k'J, 

(3.12) 



vi=l 



where 

2n 



GC"i(£;,...,*;„) -^f:-'s'"j d-.---d-2.j l^'-'Sf (n^""-] 2-4fE-< 



X Gf''^'^\k[, /c'^; /c'l, /C2n)/ai,...,a2„('^l, i^2n-l] 



X 

ai ,..,a2„=± 



(3.13) 



(i|^"^"(x2n,0)---^"n^l, 0)1 1)^=0, 

and (^(^'^+"^) ig ^Jie Fourier transform of 

(0|T ($i(x;) ■ ■ ■ ^m{.x'^)0{Xuh) ■ ■ ■ 0(x2n, t2n)) |0)g=0- (3.14) 

(For n = in ( |3.13|) the —i is omitted.) 

The integrands in ( |3.13| ) are simple to evaluate. The free Green's functions Go ai'e 
products of free field propagators for the field (f). The S-correlation functions can be 
evaluated using the algebra (|2.29| ) and (|2.48| ). As usual, dividing by Z serves to remove 



'vacuum bubbles'. We will illustrate the main features by computing the 2-point function 
to order (3^ in the next section. 
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4. Photon Self-Energy and the Dielectric Constant 

4-1 Charge-Model Computations 

Due to the overall momentum conservation, one can write 



imx.mom = n(a;,fc), (4.1) 



where 



Uiu,k) = j ^G^'Hk^k'), (4.2) 
and G*-^-* is defined in (|3.13|) . On general grounds, one expects 11 to take the form 



1 H~ o 7 o / O 7 o \ o . . . 1 . 



Ii{u,k) 



uj^ — k'^ \ ijS^ — k"^ (w^ — /c^)^ 
In quantum electrodynamics, E is called the photon self-energy. The dispersion relation is 

uo^ -k'^ - E = 0. (4.4) 

In our models, E is only a function of so the dielectric constant, defined as e{u!) = k'^/u'^, 
is given by 

= 1 - (4.5) 

The free field cp correlators appearing in ( ^.13| ) are products of propagators, which 
follow from the two point function. For the charge model, 

(0|./>(x,t)</>(0)|0),=o = / ^2_fc2lA^ + ,, ^^PH^-^)- (4-6) 

Thus, G[,'"+'^(^,^';fci,...,&2n) is a product of factors Gq^'* where 

Gl^\h,k2) = {2'.)H^^\k, + k,) I . (4.7) 

For the order contribution one needs 

G^^\k, k'; ki, k2) = G^o\k, k')G^o\ku k2)+Gi^\k, ki)G^^\k\ k2)+&^\k, k2)G^^\k\ ki). 

(4.8) 
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In the sum over a^, only ai = +, 02 = — give a non-zero contribution, and from ( |2.29|) , 
( p.48|) one has 

N 

a \S-{x2,0)S+{xi,0)\ i) = -S{x2-xi). (4.9) 

One sees that the S-correlation function here gives rise to an effective 2-point interaction of 
the photons, similar to a mass term. The first term in ( |4.8| ) gives rise to a vacuum bubble 
which is subtracted. For the other two terms in ( |4.8| ), one finds that all the integrals are 
saturated with (5-functions, and one simply finds 

n(w,A;) = ^ ,0 \o - + T^ 7^^^ (4-10) 

where 

The order /^"^ computation can be carried out explicitly. The two S-operator correlators 
that contribute to the sum are 

(i \S-{x^,0)S+{xs,0)S-{x2,0)S+{xi,0)\ i) = (^^^ d{x2i)d{x^3) (4.12) 



a \S-{X4,0)S-{XS,0)S+{X2,0)S+{XI,0)\ i) = (^-J {S{X23)S{X41) + S{X13)S{X^2)) 

N 

- 2 — 5{x23)S{^2l)S{x4l), 

(4.13) 

where d{xi2) = S{xi — X2), etc. Doing the x integrals one obtains 

n(w,A;)|/34 = — (^^J j j^j^---j^2Tcdiui + .. + uj4)Gy{k,k,ki,...,k4) 



/+_ + _(wi, W2, W3)(27r) S{ki2)S{k34) + f++ (Wl, W2, W3) 



X I — I {2ny {S{kis)S{k24) + S{k23)S{ki4)) - 4n—5{ki + /c2 + /cs + 

(4.14) 

{ki2 = ki + k2 etc.) The free Green's function one needs is 

Gi^\k, k'; fci, ^2, 4, h) = Gi^\k, k') (di^Hh, k2)G^o\k3, ^4) + 2 perm.) 

^ ^ X (4-15) 

+ [G'^Hk, ki)Gf\k', k2)G^o\k3, ^4) + 11 perm.j . 
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A diagrammatic technique can be developed for organizing the computation. One 
first draws a diagram with n + 1 unconnected hues representing G^q"'~^'^\ where each 
hne is a propagator assigned the value z/(/c^ — A'^ + ie), and one repeats this diagram 
for each possible assignment (/c, /c', /ci, /c2n) of the ends of the propagators. One then 
links them according to the 5-functions in the interactions generated by the S-correlations, 
and integrates over all remaining momenta including the factors f{u))- The S-correlation 
functions are such that they generate new kinds of interactions at each order, so we refrain 
from outlining a complete set of rules here. 

For example, the term proportional to S{ki2)S{k34) leads to the diagram in figure 1, 
which equals 



uJo(3Y fNY 1 f duj2 1 



2 J \LJ (u; + wo)(A?2 - A2)2 y 27vi {u2-uj + ie){u2-uJo + ie){u^~k^-A^ + ie) 

(4.16) 

The labels 1,2.. in figure 1 refer to uji,uj2... and the structure of the diagram implies 
ui = —u,U2 = —U4,u!s = u; this leads to the integrand ( [1.16|) . The U2 integral is easily 
done. There are poles at w — ie, ujq — ie, and ±{y/k^ + A^ — Closing the contour 

in the upper half-plane, one only picks up the pole at lj2 = — V/c^ + A^. The expression 
( [1.16|) becomes 



uoPY fN\'^ 1 f 1 



(4.17) 

There are a total of 24 terms of this type (after subtracting vacuum bubbles) in ( |4.14| ), all 



proportional to (^)^. 



w,k 1 2 4 3 -w,-k 

X X 



Figure 1. Diagrammatic representation of ( 4.16 ). 



The term in ( [1.14| ) proportional to S{ki + k2 + ks + /C4) is a new effective 4-point 
interaction of the photons. In the diagrammatic scheme described above, one contribution 
for example is shown in figure 2. This has the value 

8L^'^°^^ 2u;o(u; + u;o)(u;2-fc2-A2)2y (27r)2 (^2 _ ^2 + ,,)(^3 _ + ' > 
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There are a total of 12 terms of this kind. 

Upon summing these contributions, one finds 



where 



E4(w) 



"4 '^o 



+ 



(u;2_p_A2)3 



dudk 



2^2 



£4(0;) 



(4.19) 



1 



(271)2^2 



A;2 - A2 + ie 



2uJo 



UJ 



cuo + ^e) 



(4.20) 



The first term in ( |4.19| ) arises from the sum of aU diagrams of the kind shown in figure 1, 
whereas E4 is the sum of the diagrams of the kind shown in figure 2. Thus we have verified 
to fourth order the structure anticipated in ( [1.3|) , with E = + E2 + S4 + 0(/3^). 

In theories with more conventional perturbative expansions, the self-energy is the 
sum of one-particle irreducible 2-point diagrams. One-particle irreducibility is likely to be 
more transparent in a perturbative expansion based on the 61,2 fermion fields rather than 
directly on the 5" operators, but we don't develop this here. 

The left-over integrals in E4 identify it as a one-loop contribution. The final result 
for the self-energy of the charge-model is the following: 



E(u;) = + 



(3^ N 



2 Lu'^-ul 



X 



1 + 



Stt 



2ul 



A2 



2A^ 



log 



'ul - A2 



(4.21) 
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The size of the one-loop corrections is determined by the dimensionless parameter /3^/87r. 
The expression ( [4 .211 ) is the main result of this section. We remark that setting A to zero 
leads to infra-red divergences. 

The self-energy depends on only through the combination 13'^Nujq/L with units of 
mass^. This has a universal meaning in terms of the volume density rj of atoms: 

m = -^-jruo = - — = 27Tnuo^ri, (4.22) 

O Ij T Ij C 

with rj = limiting case of the quantum optical chain discussed above, 

~ 27ra;o P, (4.23) 

where p is density per unit length of atomic impurities. To obtain an idea of orders of 
magnitude, for hujQ = lev, one finds m = lev for p = 10^/cm. 

4-2 Current- Model Computations 

Consider next the current-model. The photon self-energy should now be defined as 
follows: 

{m<Pix.t)dt<P{Om = ^ j ^ exp(-zfc-x) ^^ -B-T. - ^^-^^^ 

The general expressions for the charge-model we obtained in this section still apply, 
except that now G\^\ki, ...^k^) is the Fourier transform of the correlation function 
{Q\dt4>{,xi) ■ ■ ■ dt4>{xn)\Q) g=o. In particular, ( [4.7|) is replaced by 

Gf\k^M) = {2n)^^^\h+k2)^^. (4.25) 

k( + ie 

Repeating the above computations, one finds 

^ , , (3^ N ujquj^ i N o f dudk uP' 
E(w) = — — —J 2 ~ ~A~ri^ ^ 



2 Luj'^-ujI 4L" V (27r)2^2_^2+^, 
2uJn 



{up' - ijjl){u) - Wo + ^e)2 



+ 



[up - u)^'^{uj - cuo + ie) 



(4.26) 

One can perform the integrals as before. An important difference from the charge- 
model is that here the one-loop integrals are ultraviolet divergent. There is a natural u.v. 
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cutoff in the model since it is implicit in ( p.4|) that the wavelengths of the photons are 
large compared to the size of the atoms. We introduce a u.v. cutoff /i as follows: 



The final result is 



2 Luo'^-ul 



1 + 



47r 



Stt 



(4.27) 



log 



(4.28) 

Here, the imaginary part arises from choosing a positive argument of the log, i.e. > f^O' 
which is physically sensible if resonant photons (with energy wq) have a longer wavelength 
than the size of atoms. Note however that Im(E) = when ij? < ujq. As we show below, 
the occurrence of this imaginary part is related to the spontaneous decay lifetime of the 
atoms. 



5. Polaritons and Spontaneous Decay 

One sees from the self-energy that the spectrum consists of two 'polariton' branches. 
To lowest order in the dispersion relation for the current- model reads 



^ {{ul + + 4m^) ± ^J{uJ^ ~ - 4m^y + 16m2u;2^ . (5.1) 



These two branches are plotted in figure 3. As m ^ 0, the two branches become cu^ = uJq 
and c<;2 = k'^^ i.e. optical phonon-like and photon-like respectively. For finite m, the co^ 
branch is phonon-like for small \k\ but photon-like for large |/c|, and visa versa for the 
LO- branch. In general one has quasi-particles with both atomic and photon degrees of 
freedom. One has uj+{k = 0) ~ c<jo + '2m'^/uJo, and uj-{k = oo) ~ ujq. Thus there is a gap 
between the two branches with 

2m2 N ^ ,2 

i?gap = = ^-=47rd\ (5.2) 

UJq 4 Li 



For the idealized quantum optical chain with fi'^/A-K ^ 1/137, if N/L = 10^/cm, then 
^gap = 2ev. 

The analagous formula for the charge-model is 

^ {{ul + P + A2) ± ^{ul -k^- A2)2 + 16m2^2^ . (5.3) 
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For A^/wq, /oJq <^ 1, the gap between the two branches is the same as in ( p.2| ). The 
main difference between the charge and current models is for the uj- branch near /c = 0. 
One finds uj'^{k = 0) ~ — Am?. For A^ = Am?, uj- = \k\ for /c 0, and the dispersion 
relations ( ^.31) and ( |5.1| ) are actually identical. In ^ the analagous problem with harmonic 
oscillator defects rather than two-level atoms is studied, and there it is shown that the 
analog of the cancelation A^ = Am^ occurs. This justifies enforcing A^ = Am^ in order to 
obtain a physical photon dispersion relation near /c = 0. 



Due to the imaginary part of E, and consequently the imaginary part of the dielectric 
constant, there is attenuation of plane waves in the medium. The origin of this imaginary 
part is the finite lifetime under spontaneous decay of single excited atomic states. To see 
this, note that for a small shift of ujq 





Figure 3. The dispersion relation u verses k for the current-model. (For = 1, 
m = .5.) 




(5.4) 
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Comparing this with ( ^4.28| ), one sees that the imaginary part of E for the current- mo del 
can be interpreted small imaginary shift to lvq: 



= I — 

ujQ 8 



(5.5) 



This leads to the decay of single atoms e^'^o* — > e^'^o^e ^Z"^, with r given in ( p.20| ). 

The attenuation constant a is defined such that the intensity of radiation decays as 



In terms of the dielectric constant, 



a 



Ime 



(5.6) 



For the current-model, one finds. 



a 



47rm — 



2\2 



8n J (w2 _ ^2) 



\UJ\. 



(5.7) 



6. Renormalization Group for the 2-Level Splitting 

We have seen that in the current-model the photon self-energy depends on an ultra- 
violet cutoff fx. A basic idea of the renormalization group is that couplings also depend on 
the scale fi in such a way that physical quantities are independent of fi. The result ( [4.28| ) 
has precisely a form that allows us to consider loq as a function of n, since 

The //-independence of E amounts to the renormalization group equation 

or equivalently 

E{u, uq{h2),H2) = S(w, uo{iJi),iJi). (6.3) 
To order Z^^, ( |6.2|) leads to the beta- function: 
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Since (|6.2| ) could be satisfied assuming no dependence of /3 on we conclude that up to 
order jS^, P is unrenormalized. 

The beta- function ( |6.4| ) means that the parameter ujq, which has engineering dimension 



(6.5) 



1/time, has an additional anomalous dimension of — /3^/47r. Integrating (^ 

^0(^2) ^ / fj^Y 



Thus, as /i increases, ujq decreases, reaching an ultraviolet fixed point at ujq = 0. 

One can derive a relation which describes the behavior of E as one scales the dimen- 
sionful parameters lj,uio,L. Ordinary dimensional analysis implies 

E(w, uJoi/i), L, n) = //^E ( -, '^^^^^ , Ln] , (6.6) 

where E is a function of dimensionless parameters. Rescaling all dimensionful parameters 
by a dimensionless parameter s, and using the renormalization group equation ( p. 3D , one 
has ^ 

^ E — , , s ] = fi S — , s] . (6.7) 

Taking = S/U, one obtains the scaling equation 

E(su;, sujq{ii), L/ s, h) = s^E(a;, uJo{sii), L, n). (6-8) 
This means that at higher energies suj, the two-level splitting uoifx) is screened to u!o{s^) = 



7. Comparison vi^ith the Reduced Maxvi^ell-Bloch Theory 

Two approximations commonly made in the quantum optics literature are the so-called 
slowly varying envelope and rotating wave approximations. In this section we compare the 
above results with the analagous results obtained in these approximations. 

6. 1 Definition of the Model 

The slowly varying envelope approximation is suitable for dealing with near resonant 
phenomena. Consider the mode expansion of the free scalar photon field: 

</>(x,t) = I {a{k)e-'^- + a\k)e^^'-) , (7.1) 



'2-K ^2\k\ 
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where k = {\k\,k). The photon creation operators satisfy 

[a{k),a\k')] =d{k-k'), 



(7.2) 



which imphes ( p.l3| ). We suppose that near resonant photons are most important, and let 

k = uq + ke, (7-3) 
where k^ denotes an 'envelope' wave vector. Letting \k\ ~ uq, one has 

(j){x, t) ^ e-^'"°(*-^V(a^, t) + e^'"°(*-"V^(a;, t), (7.4) 



where 



^(x,0) = [ 4t a(A;e)e^''=" 



V^t(^,0) = -^ [ ^a\kM 



(7.5) 



a(/ce) = a{ujQ + /Ce). (7.6) 



and 



The operators a, satisfy the same commutations as 

In classical theory, '0 and i/;^ are referred to as the slowly varying envelopes if k^ ^ Wq, 
which implies 

\d:,iP\^ojo\ipl <a;o|^|. (7.7) 



Note that in expanding the field as in ( [7.4|) , we are quantizing about a right-moving 
plane wave. One can also begin with envelopes of left-moving waves separately; however 
we will not consider interactions between the left and right moving envelopes. 

Using ( |7.7| ), one finds 



dxdt]^ {{dt(l)f - {d^(t)f) ^ 2iuG j dxdt i;\d, + dt)ip. 



(7.8) 



Without any additional interactions with atoms, the equation of motion is [dj; + dt)^^ = 
0, which means that the envelope also consists of right-moving excitations only at zero 
coupling. The canonical commutation relation which follow from ( [7.8|) are 



[^(x,t),ij\x',t)]=^5{x-x'), (7.9) 



and is compatible with (|7.5| ). 

21 



The interaction (|2.31| ) contains terms with both photon creation operators o) and 



5"^ operators which excite the atoms. Such terms lead to vacuum fluctuations wherein 
both photons and atoms are simultaneously excited, and also to real processes where 
e.g. an incoming photon excites the atom and emerges as two photons. The rotating 
wave approximation sets such processes to zero. For the charge-model, the rotating wave 
approximation, combined with the slowly varying envelope approximation leads to 



int 



(7.10) 



Note that the phases e^*'^°* in ( |7.10D cancel the time dependence ( |3.6D of 5"^ which comes 
from the hamiltonian Hq^""^. Thus we can replace e=Fit<Jo(*-a:)^± j^qt^ operators 
and set the ifg*"™ piece of the hamiltonian to zero. Thus we consider the model defined 
by the complete hamiltonian: 



H = -2iuo I dx ij^d^tjj - I dx {^S+ - tP^S') . 



(7.11) 



The first term follows from ( |7.8| ) and we relabeled the 5"^ operators back to 5"^. The 
algebra satisfied by the S"s is the same as before. In the classical context, the equations of 



motion for the model (|7.11|) is sometimes referred to as the reduced Maxwell-Bloch theory. 

The same approximations as above applied to the current-model leads to the same 
reduced Maxwell-Bloch theory ( |7.10| ), since dt{i>e~'"^°^) ^ —iuQi>e~'"^°^ . 



6.2 Exact One-Polariton States 

One can construct the one-polariton state exactly for the above model. Let us define 
momentum space S-operators as follows: 



dx 
dx 



e±'^"^±(x,0) 



e-'''''S3{x,0) 



(7.12) 



The hamiltonian now reads 



H = j dk^ ke a^{ke)a{ke) -i^^^ J dk (^a{k)S+{k) ~a\k)S- (k)'^ . (7.13) 
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The state defined in ( p.44 ), is now, in contrast to before, tlie exact ground state: 



H\Q) = 0. 



(7.14) 



By dropping ifp*"™ in ( [7.13|) , we liave merely sfiifted tlie ground state energy by Nujo/2. 
Below, we will need S3{k)\Q). From ( [TT^ ) and ( P?i5D , one has 



Ss{k)\n) 



N 



(7.15) 



where 5l is a delta function / dkdL{k) = 1, and 

5l{0) 1 



In this way. 



lim 

L— »oo L 



2n' 



|0). 



(7.16) 



(7.17) 



2 -V /I / 2 

Understanding that the polariton quasiparticle is a combination of photon and atomic 
degrees of freedom, let us take as an ansatz for the one-polariton states: 



k,) = ia^ik,) - iX{k,)S+{k,)] \n) 



where A is some function of kp. One has 



[H,aHkp)] = k, a\kp)-i(^J^S+{kp 



H, S+{k 



/ ^^(k')Ssik'-k,). 



This gives 



Thus, |/Ce) is an exact eigenstate of H, 



|0). 



(7.18) 



(7.19) 



(7.20) 



H\ke) = UJe{ke)\ke), 



(7.21) 



with 



_ f3 /I^l 

2VTa' 



(7.22) 
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if kf, satisfies 



2V 2 VA 



-A 



(7.23) 



One can thus view A as a spectral parameter, parameterizing tlie envelope energy and 
momentum We, /ce- Eliminating A one finds the dispersion relation 



m 



(7.24) 



where m? is defined in ( |4.22|) . The above dispersion relation is exact; thus one sees that 
there are no one-loop corrections of the kind computed in section 4. This also means that 
the reduced Maxwell-Bloch theory does not incorporate spontaneous emission effects. 

The result ( |7.24| ) can now be compared with the result (|4.25|) obtained to lowest 
order in perturbation theory. Recalling that and /cg denote envelope quantities, we let 
u) — uq + k = uq + ke and substitute in (|4.25| ). One obtains 



2uoUe + ul — 2uoke — k^ = Am? 



Ur 



2UeU>0 + 



Using the slowly varying envelope inequalities. 



0' 



kl<^ujl, 



(7.25) 



(7.26) 



one obtains precisely (|7.24|) . 

Remarkably, it is known that the model (|7.11 ) is integrable. The Heisenberg operator 
equations of motion are 



(7.27) 



(3uo 



(3uq 



2 ' 2 
These equations of motion have a zero-curvature representation 

[dt + Aud^ + A^] = {), 

where A^;, At are auxiliary 2x2 matrices of quantum operators: 



(7.28) 



A. 



At 



hlh2 S+ 
S- h\hi 

-i) 



+ 



1^ 

16 

1^ 
H 16 







+ 



(3UJ0 











(7.29) 
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Above, fx is an arbitrary spectral parameter, and requiring (|7.28|) to be valid for all fx is 
equivalent to ( |7.27] ). 

The zero-curvature representation allows the model to be solved by the Quantum In- 
verse Scattering Method]^, as was carried out by Rupasov[0. The integr ability leads to a 
Bethe-ansatz construction of the multi-particle states that generalizes the above construc- 
tion for the one-polariton states. 



8. Semi-Classical Analysis 

The reduced Maxwell-Bloch equations ( |7.27| ) have been extensively studied in a semi- 
classical approximation. See e.g. [0 0] and [H . The nature of this semi-classical approxi- 
mation is the following. Consider the expectation of the equations ( [7.27|) in the state 
Let us assume that the atomic and photon correlators are approximately decoupled: 

{^S+)^{i;){S+), (8.1) 

where (O) = {Q\0\Q). In this semi-classical approximation (?/') is now interpreted as a 
classical electromagnetic field. Define 

p^^S+±S-. (8.2) 

Imposing a reality condition (tp) = (tp^), one obtains the c- number equations 

dt{Ss) = -Pu;o W{p+) (8.3) 
dt{p-) = 0- 

Since {p~{t = 0)) =0, the last equation above allows us to impose {p~) = for all times. 

As a consequence of 5"^ being a Casimir for su(2), (5") ■ (5") is a constant of the motion. 
Having set {p~) to zero, this implies 



(^3)2 + = constant = (J^j . 



(8.4) 



This constraint can be parameterized by introducing an angle function 0(x,t): 



{Ss{x,t)) = cos{pe{x,t)), {p+{x,t)) = sm{(3e{x,t)). (8.5) 
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The equations (^]^) then imply 



(8.6) 



Inserting this into (|7.27|) one gets the sine-Gordon Uke equation: 



{d^ + dtd^)e 



m? N 



sin(/3e). 



(8.7) 



The classical soliton solutions to this equation were observed some time ago by McCall 
and Hahng. 

The sine-Gordon (SG) equation is easily seen to be consistent with the one-polariton 
dispersion relation obtained above in perturbation theory. Taking fi to be very small, and 
expanding the sin(/90) leads to the linear equation 



with a dispersion relation that is precisely ( |7.24|) . 

The classical SG equation has a rich spectrum of solutions consisting of solitons and 
breathers. The lowest energy breather solution can be identified with the polariton, as 
( p.8|) shows. The existence of these solutions suggests that the quantum Maxwell-Bloch 
theory may have a rich spectrum of bound states in addition to the polariton. In the next 
section we attempt to study this question by considering the quantum version of the SG 
theory. 

9. Quantum Sine-Gordon as an Effective Theory 

In the last section we saw how the classical sine-Gordon equation emerged from the 
reduced Maxwell-Bloch theory in a semi-classical approximation wherein the electromag- 
netic field was treated classically. Suppose one attempts to re-quantize the semi-classical 
treatment by quantizing the sine-Gordon theory in the canonical manner. What such 
a quantum theory has to do with the fully quantum Maxwell-Bloch theory is a delicate 
question. One can hope that an intricate manifestation of the correspondence principle in 
the end will save the day. In this section we explore these issues and conclude that the 
quantum sine- Gordon theory has some validity. 

In order to make use of the standard quantization of sine-Gordon, let us make a change 
of variables 



(a2 + aia,)e 



-m^ e. 




X = 2x — t, 



t = t, 
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such that 



+ dtd^ = ^ - 5f . (9.2) 



The action which leads to the equations of motion (|8.7|) is 



S — "f dtdx 



1 777^ 

- {d^d^ - choche) + — cos(/3e) 



(9.3) 



where 7 is an arbitrary constant. In the classical theory the constant 7 is irrelevant, i.e. 
the classical equations of motion are independent of 7. In the quantum theory however, 7 
determines the fundamental commutation relations: 

[e{xj), d^ix'J)] = ^S{x - x'), (9.4) 

and is thus meaningful. 

To promote to an operator and impose the commutation relation ( |9.4|) is potentially 
perilous given the origin of G, i.e. as a way of solving the c- number constraint ( |8.4]) . Let us 
try and interpret the quantization of G by replacing (|8.5| ), (|8.6| ) with operator equations: 



53 = cos(/3G), p+ = sin(/3G) (9.5) 

dtO = ujQ tp. (9.6) 

The S-commutation relations ( p.29|) , after setting p~ = 0, imply [Ss{x,t), p^{x',t')] = 0. 
This is consistent since [G(a;,t), Q{x',t')] = 0. 

In 1^, we argued that to connect with the quantum Maxwell-Bloch theory, one must 
take 7 = 10. Let us repeat a version of this argument here. First, note that the opera- 
tor equation (|9.6|) combined with the commutation relation ( [7.9|) does not by itself give 
a commutation relation of the kind ( [9.41 ). In a sense, one must go to next to leading 
order in the slowly varying envelope approximation in order to obtain (|9.4|) , as follows. 
Recalling the reality condition imposed in the previous section (tp) = (tp^), let us impose 
this classically. Then cj) = 2il; cos{u!o{t — x)). Since dt(j) ~ — 2t<;o'i/' sin(ci;o(i — x)), we have 
dfCp ~ —2dtOsin{uJo{t — x)). The term in the action (|2.12|) for (p that determines the 
commutation relations is 

dxdt ]^dt(t)dt(p ^4: I dxdt d^dS SYa^{uJo{t - x)). (9.7) 



^ The discrepancy of An between equation 3.4 of and equation ( p.lTD is because Heaviside- 
Lorentz units were not used throughout in Q; in particular the 47r in ( |2.8| ) was omitted in Q, 
which amounts to a redefinition of e. All formulas in this paper are in the esu system of units. 
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Averaging over the rapid oscillations let us replace sin^(a;o(i ~ x)) by 1/2. Then the 
commutation relations which follow from ( |9.7| ) are 



[eix,t),dMx\t)] = -5ix-x') 



(9.8) 



Since d{x — x') = 2d{x — x') at equal time, ( P?^ ) is precisely (|9.4| ) with 7 = 1. 

Further evidence for the relevance of the quantum SG theory defined by (|9.3| ) with 
7=1 comes from the perturbative computations in section 6, in particular the beta- 
function ( |6.4| ). From the 2-point function 



(o|e(x,t)e(o)|o) 



m=0 



--3-log(P-j2), 

471 



one has 



(0|e 



i/3e(x,t) -1/30(0) 



|0)^=o 



1 



(9.9) 



(9.10) 



(t2 - x2)/3V4^ 

This implies that the cos(/90) operator in ( |9.3| ) has anomalous mass dimension /An. 
Since the action is dimensionless, m? has a dimension of 2 — 0^ /An. Since oc ujq, 
then this precisely corresponds to the beta-function (|6.4|). 

One does not expect of course that the quantum SG theory precisely reproduces the 
quantum corrections computed in section 4 for the fully quantum Maxwell-Bloch theory. 
One can study this explicitly by computing one loop corrections in the SG model. Ex- 
panding out the cos(/3G), the lagrangian is 

2 



1 . 
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(9.11) 



The one loop contribution to the photon self-energy can be computed using standard 
perturbation theory. In the conventional coordinates x, t, the frequency and wave- vector 
cD, k are related to the envelope quantities hy ui = — kf,/2, /c = /c/2, so that ui"^ — k"^ = 
ul — Uekg. Introducing a u.v. cuttoff Jl as in section 4, one finds the dispersion relation: 

2 



ujpkp — m 



1 




0. 



(9.12) 



This should be compared to the dispersion relation cu^ — /c^ — S = for the current-model 
near resonance, i.e. with ui = uq + u^, k = uq + Ue- Using the slowly varying inequalities 



( [7.26|) , and the expression ( |4.28| ), one obtains 

32 q2 



UJeke 



m 



1 + 



An 



9>n Up \ V '^n 



ZTT 



0. 



(9.13) 
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The expressions (|9.13 ) and ( |9.12|) , including the imaginary parts, agree when both cutoffs 
are large and when Ue ^ ujq; however this contradicts the slowly varying inequalities. 

The conclusion of the above analysis is that the quantum SG theory captures some 
aspects of the fully quantum Maxwell-Bloch theory, in particular the current-model defined 
in section 2, and does incorporate spontaneous emission, but is not equivalent to it even 
in the slowly varying envelope approximation. 

As a step toward understanding the spectrum of the quantum Maxwell-Bloch theory, 
one can assume the approximate validity of the quantum SG description. The quantum 
SG spectrum is known to consist of breathers, the lowest mass breather being the particle 
associated with the SG field O itself, and a pair of solitons. See e.g. The mass of the 
n-th breather is given by 

™„ = 2™,si„(^j, n= 1,2, ...<-, 

where is the mass of the soliton. As /3 ^ 0, the mass of the lowest breather mi 
approaches m thus the polariton is identified as the lowest breather. The higher breathers 
are polariton bound states. As /9 — > 0, TJig ~ 8m//9^, thus for very small the mass of 
the soliton can be very large compared with the polaritons. In the quantum theory, the 
polariton can actually be viewed as a bound state of two solitons. 

It would be very interesting to understand whether aspects of this spectrum can be 
seen in the models defined in section 2. 



10. Conclusions 

We have defined some models which describe quantized radiation in interaction with 
a medium of two-level atoms arranged in a fiber geometry. Our main computational 
results are the photon self-energy ( [4.21| ) and ( [4.28| ), which determine the first quantum 



corrections to the polariton dispersion relation. We also compared our results with known 
semi-classical results in the slowly- varying envelope and rotating- wave approximations, and 
argued for the approximate validity of the quantum sine-Gordon theory. We found that 
the model which follows from an interaction hamiltonian —d ■ E (current-model) is better 
behaved than the model which follows from the minimal coupling p ^ p — eA, the latter 
suffering from infrared divergences. In the current model we derived a renormalization 
group equation for the energy splitting of the two-level atoms which follows from the beta 
function (|6.4|), and shows that the splitting is screened at higher energies. 
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Though the quantum corrections are generically small, we hope that the trend toward 
fabricating smaller optical devices will eventually lead to the observation of these quantum 
effects. 

The models defined in section 2 deserve further theoretical study, in particular it would 
be interesting to determine whether they have a bound state spectrum that resembles the 
spectrum of the quantum sine-Gordon theory. 
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